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The influence of the nuclear medium upon the internal structure of a com- 
posite nucleon is examined. The interaction with the medium is assumed to 



depend on the relative distances between the quarks in the nucleon consistent 



with the notion of color neutrality, and to be proportional to the nucleon den- 
sity. In the resulting description the nucleon in matter is a superposition of 
the ground state (free nucleon) and radial excitations. The effects of the nu- 
clear medium on the electromagnetic and weak nucleon form factors, and the 
nucleon structure function are computed using a light-front constituent quark 
model. Further experimental consequences are examined by considering the 
electromagnetic nuclear response functions. The effects of color neutrality 
supply small but significant corrections to predictions of observables. 
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I. INTRODUCTION 



Nucleons are composite color singlet systems made of quarks and gluons. Their wave 
functions consist of many configurations as illustrated in Fig.|l]. Some configurations are 
simple with only three current quarks, most are more complex with many partons. Configu- 
rations in which three quarks are close together have been dubbed point-like configurations 
(PLC). PJj According to perturbative quantum chromodynamics (pQCD), such configura- 
tions are responsible for high momentum transfer elastic scattering reactions . Further- 
more, the effects of gluon emission from closely separated color-singlet systems of quarks and 
gluons tend to cancel. This means that for processes in which one adds amplitudes before 
squaring, i.e. coherent processes, point like configurations do not interact with surrounding 
media. 

That color neutrality can suppress interactions is similar to the cancellation of inter- 
actions which would occur if an electron and positron would move together at the same 
position through a charged medium (see for example Ref. Q). The color neutrality feature 
of QCD has been verified. It is responsible for the scaling of structure functions at low, but 
not too low, values of the Bjorken scaling variable xbj- (See the reviews 0-0].) Furthermore 
hadron-proton total cross sections 07^ grow linearly with the mean square radius of the 
hadron [|8]]. 

There are also configurations in the nucleon wave function in which the partons occupy 
a larger than average size. We call these configurations blob-like configurations BLC or 
huskyons |§. These configurations have complicated strong interactions with the medium. 
(See [[| which explores the consequences of huskyons.) Here we model the interaction of 
the nucleon with the surrounding medium in terms of the inter-quark separation within 
the nucleon to investigate the role of color neutrality and its consequences in nuclear wave 
functions. 

Sixty years of studies of nuclear properties have shown that the nuclear wave function is 
dominated by clusters of color singlet objects with the quantum numbers of nucleons and 
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mesons. The success of the nuclear shell model is a testimony to this fact. It is certainly 
possible, however, that although bound nucleons largely maintain their identity, bound and 
free nucleons are not identical. Indeed, there are a number of interesting experimental 
findings which may indicate that medium modifications of nucleon properties are relevant. 



One of the most spectacular examples is the observation of the first EMC effect JTD 
showing that the structure function of the bound nucleon was suppressed at large XBy 
There have also been numerous studies of the (e, e') reaction which find that the longitudinal 



response function Rl is suppressed, while the transverse response Rr is not. See Ref. [IT 



and references therein. The suppression of Rl is natural in some theories [12| . This lore has 



been challenged recently by Jourdan |]T3| who argues that the "so-called quenching is mostly 
due to the limited significance of the data" and that including data at high energy loss 
uj leads to the result that "no A-dependent quenching is observed". However, the various 
errors listed in Jourdan's Table 2 allow for up to about 15% effects. 

It has also long been stated that the nuclear value of the axial coupling constant Ga was 



less than its value in free space WM . Furthermore the large pionic enhancement expected in 



many models of nuclear structure was not observed [TJ| , but there may be hints of a pionic 



enhancement IITB 



Thus there are several indications that nucleons could be modified by the medium. Since 
color is at the heart of QCD, one might suspect that looking at the consequences that color 
neutrality might have for bound composite nucleons might be worthwhile. These are the 
presumably small effects of nucleonic polarization. Nevertheless, as we show here, these 
effects can have nontrivial consequences. Before describing in detail our treatment of the 
nucleon in the medium, it is necessary to mention some of the possible objections and also 
to discuss some of the relevant history. 

The most significant objection is that there could be many effects other than color neu- 
trality which cause significant modifications of the properties of nucleons. This is certainly 
a valid point. For example, in the Walecka model ]I7J the medium provides a decrease of 



the nucleon mass by as much as 30%. See also the quark-meson coupling model of Ref. |T8 



and the works of Ref. |T9| . |20|| . It is entirely conceivable that similar effects to the nucleon 



form factors could occur from a purely hadronic mechanism, i.e., independent of the relative 
quark separation. It has also been argued that a more general scaling of hadronic masses 



may arise in the nuclear medium. [21 1 The extension of QCD sum rules to finite density 



further confirms the presence of large scalar and vector self-energy corrections P2[ , P3"| . Nev- 
ertheless, these effects have a different origin than that of color neutrality, and our aim 
here is to simply study the consequences of this single effect in several different situations 
by identifying a pattern of predictions. The study of many reactions is necessary, so we 
derive a formalism applicable to many reactions. This formalism may be useful in the study 
of other mechanisms for medium modifications as well as other corrections to the impulse 
approximation. 

Another objection is that it is not clear when pQCD, and hence its discussion of point-like 
configurations is applicable. The question of how high the momentum transfer must be for 
pQCD considerations to dominate calculations of form factors is controversial. The work of 
Refs. is merely the first example of a long history. More data are needed to settle this 

question. However, point-like configurations may arise from non-perturbative considerations 
as well PB|P?|. The singular nature of the non-perturbative confining interaction can lead 



to point-like configurations. Indeed, many constituent quark and Skyrmion models seem to 
indicate the presence of point-like configurations. It is therefore of interest to consider an 
interaction with the nuclear medium which accomodates such objects. 

Another possible objection is the failure to observe clear and convincing evidence for 
color transparency CT. This is the idea that a high Q 2 quasielastic reaction produces a 
point-like configuration which does not interact with the surrounding nuclear medium. In 
this way initial and/or final state interactions are reduced. A BNL experiment ]2B| shows 
some evidence for CT in a (p,pp) experiment at beam momenta of 6, 10 and 12 GeV/c (4.5 
GeV 2 < Q 2 < 8 GeV 2 ) but the (e, e'p) NE18 (Q 2 = 1, 3, 5, 7 GeV 2 ) experiment at SLAC |2| 
sees no such evidence. A recent FNL exclusive p production experiment saw evidence for 
color transparency [30|. The Q 2 was up to about 9 GeV 2 . The BNL and FNL experiments 



may not have guaranteed the necessary quasielastic nature of the experiment. The SLAC 
experiment did that. One interpretation of these results is that point-like configurations 
are produced but expand before leaving the nucleus. This expansion is a bigger effect for 
the NE18 experiment, which has the lowest energy outgoing protons. Expansion is also a 
big effect for the BNL experiment, but the incident proton has high enough momentum 
for this expansion not to completely kill the influence of the point-like configurations. The 
expansion effects are smallest for the FNL experiment, but the experimental resolution may 
be a problem. 

None of the considerations of the above paragraphs are sufficient to rule out the existence 
of point-like configurations, or the effects of color neutrality in nuclear physics. Indeed, 
point-like configurations were introduced in the context of the EMC effect |l] . Frankfurt and 
Strikman postulated that such configurations do not feel the attractive nuclear interaction 
and therefore appear in a smaller percentage in nuclei than in free space. This leads to a 
reduction of the nuclear structure function at xbj ~ 0.4 — 0.6. This effect also reduces the 
cross sections predicted for (e, e'p) reactions, making color transparency harder to observe 
[5I|. Frankfurt and Strikman treat their effect by taking the suppression to vary essentially 



as a theta function in Q 2 , i.e., present at Q 2 above a certain value but completely absent 
for lower values. Desplanques finds a similar term by including the spatial variation of the 
nuclear scalar and vector mean-fields over the volume of the nucleon |[32|| . Desplanques's 
treatment was non-relativistic. Note also that Jain and Ralston |33| have argued that hard 
processes are modified significantly by medium effects. 

Our aim here is to present a method for doing calculations which may be applied both 
at high and low values of Q 2 . We also obtain a Lorentz covariant formulation. Here is 
an outline. Our model and basic formalism are presented in the next section. Section 3 
deals with some simple examples which allow us to study our approximations and expose 
the need for a relativistic formalism. The relativistic formalism, which applies the recent 



work of Schlumpf [|54)-|3"S]|, is discussed next. Section 5 is concerned with a presentation 



of our detailed numerical results. We study medium modifications of the electromagnetic 



proton form factors Ge and Gm and the axial vector form factor Ga relevant in the weak 
interaction. We discuss consequences in the e, e' and e, e'p reactions. A concluding section 
discusses and asseses these results. 

II. FORMALISM AND MODEL 

We begin from general considerations. Suppose a color singlet baryon moves in the 
nucleus subject to a Hamiltonian H given by 

H = H + H 1 (1) 

where H\ is a perturbing Hamiltonian, and 

H = H^ m + H^ 1 . (2) 

Here Hq™ and H^ el describe respectively the motion of the center of mass (including effects 
of the medium on the center-of-mass motion), and the relative motion of the internal degrees 
of freedom independent of the medium and the center-of-mass motion. The perturbation Hi 
describes mixing between the center-of-mass and relative motion, and thereby incorporates 
the effects of the medium on the internal wave function. 
We shall assume that H 1 is separable, that is, 

Ht = Hl m h\ e \ (3) 

where iff 71 and h™ 1 act only in the center-of-mass and relative motion spaces respectively. 
Our separation of the Hamiltonian in (|l|) and (0) entails that h{ el requires a change in the 
internal structure of the baryon. This implies that h[ el has no diagonal elements, and that 
the baryon in the medium is described as a superposition of excitations. 

We wish to consider the eigenstates |\l/ nm > of the hamiltonian if in ([]]), defined by 

11 | ^ nm J -'nm \ * nm ^ ■ V / 
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The form of H Q in (|2|) allows the unperturbed states to be written as a direct product of 
center-of-mass and relative motion states as 

\* ( ± >= |*» > \<t>m > ■ (5) 

Here > satisfies the eigenvalue equation 

where # cm |$„ >= £ n |$ n > and H r Q el \(j) m >= e m \(j) m >, with £<W = £ n + €m . 

The state vector of the full Hamiltonian can be written to first order in Hi as 

Itf >= |*(°) > + V - *ki\Hi\V<£l> ,y(fi) (J) 

\^nm T nm ^ ~ ,-,(0) „(0) \ kl ^ • \' J 

To simplify our notation we restrict our discussion to the nuclear ground state. For this case 
Eq.([?p can be rewritten as 



1*00 >= w > +E I** >< **iJ*ri*o > E J; 1 Z(o) fa > ■ ^ 

k 1^0 -&00 - &kl 



It should be noted that in (|8|), since h™ 1 has no diagonal elements, I 7^ and the sum on 
k is unrestricted. The energy denominator in Eq.([8|) is dominated by nucleonic excitation 
energies. These, which are typically hundreds of MeV, are much larger than the nuclear 
energy differences which are typically tens of MeV. Thus we write 

Ei$-E$>*e -e, (9) 

This allows the sum on k to be performed using completeness so that 



*oo >= 



1^0 e o ~ e i 



|$ > • (10) 



We define the quantity in brackets as \<p > 



10 100 > +E < mel[h > i^ > Hr, (ii) 

i#o e ° - e i 

which can be regarded as the modified nucleonic wave function. Note that this wave function 
depends on the coordinates of the entire nucleus through the operator H^ m . The use of 



Eq.fllTD allows one to evaluate the effects of color neutrality for finite nuclei; this equation 
is our principle formal result. 



The evaluation of Eq.(|TT|) depends on knowing the wave functions corresponding to the 
Hamiltonian H™ 1 . At the present time there is no complete relativistic treatment available, 
but progress has recently been made in this direction [R^J. We wish to obtain an alternate 
method of evaluating the influence of the nuclear medium on nucleonic wave functions. Such 
can be obtained using a closure approximation in which all of the strength is assumed (on 
average) to lie at an average excitation energy. We use eo — ei = AE. The expectation value 
of an operator O can then be written to first order in Hi as 

TTcm 

< 0|0|0 >=< O |0|0 O > +^ < 0o I {h[ el , O) |0o >, (12) 

where {A, B} = AB + BA. The right side of (0) is independent of the excited state wave- 
function, and can therefore be evaluated based on knowledge of the ground state wavefunc- 
tion and the excitation energy. Thus the uncertainty in the knowledge of the wave functions 
is replaced by the uncertainty in a single parameter, AE, and by the further assumption that 
AE does not depend very strongly on the operator O, i.e. that AE is essentially process 
independent. The accuracy of this approximation and assumption is investigated in section 
4. We shall discuss reasonable values of AE after discussing our choice of the operator h"[ el . 

Equation ( fLlp has a simple physical interpretation. Consider, for example, the case where 
the operator O is the electromagnetic(EM) current j M associated with a constituent of the 
state |0o >■ A modification of the form factor will occur by the mechanism illustrated in 
FigfJ. The current acting on the ground state excites an internal mode. The excitation 
decays back to the ground state via the interaction with the medium provided by h[ el '. The 
extent to which the form factor is modified is determined by the ability of the EM current 
to produce an excitation and the ability of the medium to absorb its decay and return 
the composite particle to its ground state. The modification is directly accountable to the 
density of the medium in the vicinity of the interaction through the potential, Hl m (R), and 
is suppressed by the excitation energy AE. 
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Medium effects are often estimated by taking the nucleus to be infinite nuclear matter, 
or by using a local density approximation. Eq.(^) shows that instead one can evaluate the 
appropriate nuclear matrix which depends on iJf™. 



III. COLOR SCREENING 

Consider, now the motion of a composite color singlet baryon through a nucleus. We 
are here interested in the properties of the ground and low lying nuclear states, so that the 
baryon is part of a bound state wave function. The configurations of the baryon are pictured 
in Fig.|I]. Let the displacement of the center of mass of the baryon from the nuclear center 
be denoted as R. The interactions between such a complicated system and the remainder 
of the nucleus must depend on the positions fj of the partons inside the baryon. 

In general the interactions between a nucleon and the rest of the nucleus are compli- 
cated. Here we are studying the presumably small effects of nucleonic polarization, so we 
concentrate on the necessary modifications of the central part of the nuclear shell model 
potential. This is the largest interaction to consider, and a small modification of it might 
have significant consequences. Thus we have the central potential V = V(fi,R). What else 
do we know? Color neutrality tells us that V = when fi = R for all partons i. Further- 
more interactions vanish as r 2 , where r 2 = £j<j(r*j — fj) 2 p8|-^0|. Finally we note that the 



nucleonic average over f should correspond to the standard nuclear shell model potential. 
These considerations allow us to write 

2 

V(r,R) = Vop(R)^^, (13) 

where < r 2 >=< <f)o\r 2 \<j) > is the nucleonic expectation value of the operator r 2 , Vq ~ 
— 50MeV and p(R) is the nuclear density normalized so that p(R = 0) = 1. 

We are concerned with fluctuations, so that it is convenient to rewrite the central shell 
model potential as 

2 2 

V(r,R)=V p(R)+V p(R) T > , (14) 
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< r 2 > 



so that the operator Hi of previous sections can be identified as 



r 2 — < r 2 > 

Hl = V 0P (R) <r2> . (15) 

This is the simplest form of Hi that we can write, which is consistent with known properties. 
The behavior for large values of r 2 is simply a guess, but it is reasonable to expect that large, 
blob-like configurations should have strong interactions with the medium. 

The equation includes the effects of color neutrality on the total nuclear mean field 
and makes no distinction between the scalar and vector mean fields. The effects of color 
neutrality on the separate fields can be examined using the quark-meson coupling model 
P~8|] , but this would require the extension of the applicability of that model to high values 



of the transferred momenta. 

We note that Hi has the interesting property that 

<0o|#i|0o>=O, (16) 

which means that the color neutrality effect governed by our perturbing Hamiltonian does 
not give a first-order shift in the mass of the baryon. Such effects are by definition contained 
in the first term of Eq.(p^). 

Let us discuss some of the implications before presenting the evaluation of specific models. 
Firstly, consider the value of AE. We see that the operator r 2 — < r 2 > excites the breathing 
mode of the nucleon. Thus it is reasonable to associate AE with the the energy at which 
the first resonance of the nucleon occurs. This is the Roper resonance with AE=-500 MeV. 
We are concerned with how certain matrix elements are influenced by the presence of the 
nuclear medium, so it is convenient to define 5 < O >=< 4>\0\4> > — < 4> \O\(j)Q >. Then 



using Eq.([T5|) in Eq. (|TJ) leads to 



This shows that the nucleon's properties as measured by 5 < O > depend on the position of 
the nucleon. The factor 2Vop(0)/A£' is about +0.2, which shows that the effects we study 
here are small but not insignificant. 
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It is useful to consider examples. First suppose O = r 2 . Then 

2 „ V p(R) ± . r 4 - < r 2 > 2 , 1 
*<r 2 >=2 J ^<U <r2> |0o >■ (18) 

Writing r 4 = r 2 J2n \4>n >< (fin\r 2 , subtracting the ground state contribution, and examining 
the remainder and using completeness shows that 5 < r 2 > > 0. Color neutrality leads 
to an increase in the mean square radius of the nucleon. Note that the same argument can 
not be used to show that 5 < r 4 > > 0; indeed this quantity could be negative. Color 
neutrality does not correspond to a general scaling of the nucleon wave function. Each 
matrix element must be worked out independently. 

Another example to consider occurs in the non-relativistic quark model if O — jj, = 
52i=i,3 /^(O- I n this case the expectation value of O for a spin up nucleon is the magnetic 
moment of free nucleon. In the non-relativistic quark model the wave function is a product 
of space and spin-isospin functions. In this case, the integral over the spatial coordinates 
vanishes and 5 < fi >= 0. Similarly, in the non-relativistic quark model the modification 
to the electric and magnetic form factors must be the same. It is necessary to extend these 
considerations to relativistic models. 

It is very interesting to consider the case when O = O = S[{fi — f 2 ) / y2\S[(fi + r 2 — 
2^*3) / Vo] for which all of the quarks are at the same position. The expectation value of Oq is 
the square of the wave function at the origin. In this case the assumption that V(r = 0, R) 
vanishes is sufficient to give a result, i.e., no additional assumption about the dependence 
on r 2 for large r 2 is needed. One finds immediately (using Eq.(|17D for example) that 



5<O >=- 2 -^\M0)\\ (19) 
which corresponds to about a 20% reduction of the square of the wave function at the origin. 



This result is in good agreement with a recent QCD sum rule calculation by Jin et al. [II 

It is worthwhile to compare our approach with that of Frankfurt and Strikman Jl| . Those 
authors consider the nucleon \N > as a sum of various configurations, \N >= \PLC > + • • •. 
At high Q 2 the form factors are assumed to be dominated by the \PLC > component. In 
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the medium \PLC > is replaced by (1 + j^)\PLC >. However, r 2 \PLC >= 0, so that the 
PLC acquires a factor of (1 — Vo £ffi ) ■ This leads to 10 — 20% reductions in form factors, if Q 2 
is greater than some value large enough for \PLC > to dominate. This result is recovered in 
our approach if we assume that the only dependence of the form factors on the (modification 
of the) nucleon wave function is through the wave function at the origin. Even more striking 
is the feature that the \PLC > is suppressed, so that one expects valence quarks to carry less 
momentum in the nucleus. This leads to an explanation of the EMC effect. Our approach 
is motivated by the original work of Frankfurt and Strikman. However, we do not make the 
"all or nothing " assumption in which \PLC > suppression turns on abruptly. Thus we are 
concerned with making more detailed evaluations than and testing the assumptions of the 
early work. In particular, we wish to learn if these suppression effects persist if one uses 
more detailed models. 



IV. TOY MODELS 

To gain some insight into the results obtained for the nucleon in the relativistic frame- 
work to follow in section 5, it is useful to apply the formalism of the previous two sections to 
some simple non-relativistic constituent quark models involving (mainly) only two quarks. 
These simplifications enable us to obtain exact solutions and therefore to assess the validity 
of our approximation scheme. We shall start by using the harmonic oscillator model. How- 
ever, such a model is not expected to be a reasonable guess for large momentum transfers. 
Indeed previous work [^,0 showed that the harmonic oscillator does not have a point-like 
configuration. Therefore we also consider the case of two quarks bound by a Coulomb po- 
tential. This model assumes that the attractive nature of the color electric force dominates 
over the longer range confining force. Finally we consider a model with both a harmonic 
confinement and Coulomb term as in p6| , |27[| . 
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A. Harmonic Oscillator 



Here we start by considering the "nucleon" to be made of two quarks, and employ the 
nonrelativistic constituent quark model with H™ 1 given by the harmonic oscillator form 

H r el = f£ + \^y. (20) 
We study the nuclear medium- modified nucleon wave function <j> at R — 0. This is given by 



+ -nujy + v oP (o) 



E<j>. (21) 



2/i ' 2 r " u " ' ' ur ^-'< r 2 >_ 
We may easily obtain the exact solution by realizing that in the medium the frequency u 
is replaced by uj with 

2y oP (o) 

LO 9 

fi < r 2 > 

= ul + ^op(0)cu . (22) 

For typical values of AE = -2lu ~ -500 MeV and V = -50 MeV, we find that uj = 214 
MeV and the expectation value of r 2 increases by a ratio of 1.17. So, at the nuclear center, 
the effect of color neutrality can be significant. The expectation value of r" 1 varies as the 
inverse of the square root of the frequency, so that its expectation value is decreased by a 
factor of 0.93. 

We may also compute the form factor for this model. This is the Fourier transform of 
the square of the wave function as a function of the momentum transfer Q = \Q\ divided by 
two: 

F (Q 2 ) =< o |e i( 5 %o > (23) 
F(Q 2 )=<0|e^|0>. (24) 

Q 2 Q 2 

The free form factor F (Q) = e 16 ^o 5 while the medium modified form factor F(Q) = e . 
We see that 

Mm = 0, (25) 
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so that huge effects at high Q 2 are possible. However, this model is not realistic at high Q 2 . 

We may also consider the effects of color neutrality for very dense nuclear systems by 
increasing p(0) above its value of unity in normal nuclear matter. Eq.(|22|) tells us that uj 



vanishes for densities about four times nuclear matter density for which p(0) = |p£ = 3.75. 
Thus there is a "deconfinement" phase transition inherent in our model. However, the model 
is built on the assumption that the shell model, with its non-overlapping nucleons, is a valid 
starting point and therefore should not be applied to the situation for which the interparticle 
spacing is less than the diameter of a nucleon. 

We note that the extension of this model to the 3-quark system yields essentially the 
same results. This is because the harmonic interaction r 2 = J2i<j(u — rj) 2 = 3(p 2 + A 2 ) 
where p = (r*i — r^^l and A = (r 1 + r 2 — 2r* 3 )/v / 6- Using the same operator in H 1 once 
again leads to the result that for harmonic oscillator models the frequency in the nuclear 
medium is less than that in free space. 

Let us examine the accuracy of our first order (fo) treatment of Eq.(|I7]). In first or- 
der S < r 2 >/ = ~ 2 ^,^°" > < r 2 > compared with a model exact result of 5 < r 2 >= 
((1 + ^f^ " 1 ) -1 ^ 2 — !))• With our typical parameter, the first order increase in the mean 
square radius is 13%, while the exact result is a shift of 17%. 

The EM form factor in the medium obtained from Eq.(|l^) is given by 

V Q 2 



1 - 2- 



F (Q 2 ), (26) 



AE 6fiu 

where F (Q 2 ) = exp(—Q 2 /\6iiuJo) is the free form factor. Eq.(|26|) compares favorably 
with the exact solution with corrections of order (Vq/AE) 2 . Some caution should be ex- 
ercised, however, in the application of the approximate solution to large momenta, i.e., 
Q 2 ~ 6fiu AE/Vo, where the coefficient of the (V /AE) 2 correction becomes significant. It 
is also important to note that the interaction with the medium has produced an increase in 
the mean square radius, < r 2 >^ Tq[1 + ^(Vo/ AE)], as is reflected in both the exact and 
the approximate form factor. 
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B. Coulomb binding 



The harmonic oscillator is not expected to describe situations involving high momentum 
transfer. Indeed we may recall the asymptotic expansion: 



lim F (Q 2 



^7T d(f>l 



Q 4 dr 



+ 



16vr d 3 



r=0 



Q6 dr*™ 



(27) 



r=0 



This power-law dependence gives larger results than the Gaussian form and holds unless the 
potential is an analytic function of r 2 . The only relevant example of such is the harmonic 
oscillator force of the previous sub-section. It is useful to obtain the general expression for 
the medium-induced change in the form factor AF(Q 2 ) = F(Q 2 ) — F (Q 2 ). To first-order 
in Hi this is given by 

< O |e^'%„ >< 0n|#l|0O > 



One may then define the state vector \\ > such that 

(E - H r d )\ X >= (1 - |0o >< folWfo >, 



(28) 



(29) 



with < x|0 o >= 0. Then 



AF{Q 2 



2 < <p \e lQ -i\x >, 



(30) 



and 



167T d 3 

„ + Q 5 "* 3 



(31) 



r=0 



Our purpose here is to compare the exact first-order result of Eq. ( P9|) with that of the closure 
approximation of Eq . (|12"D . 

The above equations are general, but we can gain some understanding if we specify to 
the Coulomb Hamiltonian 



rel _ P_ 
2/2 



- e 2 /r, 



(32) 



15 



with e 2 = 4a s /3 to simulate the strong color electric force. The ground state wave function 
is given by 

, , s 1 -JL 



Mr 



—e a o . 



(33) 



with ao = l/(e 2 /i). The form factor for this state is 



F (Q 2 



1 + 



Q 2 a 2 
16 



(34) 



The Hamiltonian of Eq. (|32|) can be used in Eq . fl29|) to obtain a solvable differential 
equation. The result is the function x(r) given by 

V p(0)nal /2 



X(r) 



-e-^[ll/2-(r/a ) 2 -(r/a ) 3 /3]. 



(35) 



One may use Eqs.(|31|) and (|35|) to immediately find that, at large values of Q 2 , AF{Q 2 ) 
varies as Q~ 4 and is negative. In this model, the high Q 2 form factor is dominated by the 
point like configuration, so we may say that the point like configuration is suppressed in the 
nuclear medium. It is not difficult to use Eqs. ( pOf ) and ([35]) to obtain an exact expression 
for AF(Q 2 ). However it is more instructive to use the asymptotic expansion of Eq.flTFp to 
get AF(Q 2 ) as an expansion in powers of Q~ 2 . Then one finds 



lim AF(Q 2 

Q 2 — >oo 



V oP (0)n 



544 



(36) 



3 a 2 , Q 4 

Recall that Vq is a negative quantity, so that AF is also negative. The closure approximation 
to the change of the form factor, AF c i os (Q 2 ) can be obtained from Eq . (|i~2D as 

2V»p(0) 



which may be evaluated asymptotically as 

AF cU Q 2 - 8 ^ (0) 



-4 V ^)F (Q^ 
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(37) 



(38) 



3AE \Q 2 al) 

It is clear that the closure approximation can not be valid unless AE is taken to be a function 
of Q 2 . We may equate the AF' s of Eqs.(|5BD and ([J6|) to determine the "correct" value of 
AE. The asymptotic result is 
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A E = -^p[llQ 2 a 2 -Q8}- 1 . (39) 

This means that the magnitude of the AE decreases as Q 2 increases, so that using the 
closure approximation with a fixed value of AE underestimates the effects of the medium 
modifications. This is the principal result of this subsection. 

C. Harmonic Oscillator plus Coulomb 

A more interesting model is one which includes both a confining and a Coulomb type 
1/r term in the Hamiltonian. Thus we consider 

< - i + \^ - It (40) 

This model was considered in Refs. [p6|,p7| in which it was shown that including the Coulomb 



term leads to PLC dominance of the form factor. This occurs even though the Coulomb 
term causes only a small effect in the computed energy. 

For the calculations performed here the parameters are taken to be a s = 0.1, \i = 
300MeV/c 2 and Huq = 390MeV. Fig.[| shows how the ground state wave function in free 
space is influenced by the 1/r term in the Hamiltonian. The attraction enhances the wave 
function at the origin and changes the shape away from the Gaussian. This causes a 1/Q 4 
behavior in the form factor. Next, Fig. f| compares the full wave function with and without 
the effects of the medium. We see that the medium causes a reduction of the short distance 
wave function. The closure approximation to the wave function in the medium is compared 
with the exact calculation in Fig. [5]. The numerical results for the medium-modified form 
factor are shown in Fig.^ along with the closure result and the free form factor, from which 
one concludes that closure works well and in fact slightly underestimates the exact result. 
This is consistent with the conclusions of the previous subsection. 

It is also useful to try to understand these results using analytic techniques. This can be 
done if one works to first order in a s . First consider the free case. One may find the form 
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factor using an equation like Eq. (p8|), but with the operator Hi replaced by — 4a s /3r. The 
result is 



with 



where 



(41) 



/(*0 = E - 



i i i 



j n 2n + 1 n! 



(42) 



= Q 2 /(16(JUjj ). 



(43) 



The second term of Eq.(|4T|) leads to a (hujq) 2 /Q a behavior of the form factor which dominates 
over the exponential term for values of z Q greater than two or so. Then the free form factor 
varies as u>^ 2 . The medium modification is quite easy to implement. Simply replace u by 
u of Eq.(0)in the formulae (gl])-©. We have seen in Eq.@ that u = 0.8Qu so that 
the form factor is reduced by a factor of 0.79. This is a significant suppression. The results 
of this perturbative analysis are shown in Fig.|7|. Comparison with Fig.^ shows that the 
perturbative analysis is in good agreement with the full result. 



V. RELATIVISTIC CONSTITUENT QUARK MODEL OF THE NUCLEON 



The influence of medium modifications on computed matrix elements and form factors are 
examined using non-relativistic models in the previous sections. The need for a relativistic 
formulation is clear; we wish to compute form factors at momentum transfers Q 2 > 1 GeV 2 , 
and we wish to be able to distinguish between electric and magnetic effects. 

Here we study the effects of color neutrality on the light front nucleonic model wave 
functions of Schlumpf |]3"4|-|36|j . The first formulation of such a light front relativistic quark 



model was presented by Terent'ev and Beretskii [ f4^ . f43"f . Many authors p4]- |52| . |37| , |53| -p^, ?] 
have contributed to the development of this model. We use Schlumpf 's model because his 
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power-law wave functions lead to a reasonably good description of the proton electromagnetic 
from factors, Ge and Gm, at all of the Q 2 where data are available f34|^36| . 

In a quantum mechanical relativistic theory the commutation relations between the ten 
generators of the Poincare group must be respected. The light front approach is distinguished 
by the feature that the maximal number of seven generators are of kinematical character 
(do not contain interactions). Another feature involves the use of the light front variables 
p± = p° -|- p 3 ) go that the Einstein mass relation p^ = m 2 can be expressed as 

P~ = ((P± + m 2 )/p + , (44) 
where p± = (p l ,p 2 )- Susskind |58] noted that this equation is similar to the non-relativistic 



kinetic energy if one interprets the variable p + to be a relativistic version of the mass. Thus 
one can use relative momenta for systems involving several particles, with the result that 
the wave function is a simple product of a function involving only relative momenta with 
a separate function carrying information about the motion of the center of mass. One may 
also employ the Melosh transformation [5I| to construct states that are eigenfunctions of 



the total angular momentum and its third component. 

The light front dynamics have another important feature, stressed in Refs. PJ4"1| , that 
the diagrams with quarks created out of or annihilated into the vacuum do not contribute. 
Furthermore, one need only consider three quark components of the nucleon, if one is com- 
puting the matrix elements of "good " operators . 



Schlumpf's model is well-documented p4]-p6|. We reproduce the relevant features here 
for the sake of clarity, following his thesis closely. It is necessary to express the ten gen- 
erators of the Poincare group P M and in terms of dynamical variables to specify the 
dynamics of a many-particle system. The kinematic subgroup is the set of generators that 



are independent of the interaction. There are five ways to choose these subgroups [|6T 
Usually a physical state is defined at fixed x°, and the corresponding hypersurface is left 
invariant under the kinematic subgroup. 

The light-front formalism is specified by the invariant hypersurface x + = x° + x 3 = 
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constant. The following notation is used: The four-vector is given by x = (x + ,x~,x±), 
where x^ = x° ± x 3 and x± = (x 1 ,^ 2 ). Light-front vectors are denoted by an arrow 
x = (x + ,x±), and they are covariant under kinematic Lorentz transformations ||62|| . The 
three momenta pi of the quarks can be transformed to the total and relative momenta to 
facilitate the separation of the center of mass motion [F3] as 



B , Pi pf+pi 

P = Pl+P2+P3, ^=—\—, + > V 



Pl+Pl ' P + 



QX = (1 - £)P1± - &2± , K± = (l-T)) (piX + P2±) - VP3± 



(45) 



Note that the four-vectors are not conserved, i.e. pi + P2 + P3 7^ P- In the light-front 
dynamics the Hamiltonian takes the form 

H = ^±^, (46) 
2 p+ ' v > 

where M is the mass operator with the interaction term W 

M M ■ If . 

M ^^l_ + Ml + J±, (47) 

rj(l — 7]) 7] 1 — 7] 



2 = 01 , m l , m 2 



with mj being the masses of the constituent quarks. To get a clearer picture of M we 
transform to (73 and K 3 by 

£ 2 + g 3 #12 + K 3 



Ei + E 2 ' E\2 + S3 

(48) 



^i/ 2 = (q 2 + m 2 /2 ) 1 / 2 , E 3 = (K 2 + m 2 ) 1 /2, E 12 = (K 2 + M 3 2 ) 1 



/2 

j 



where q = (qi, q 2 , q 3 ), and K = (iTi, K 2 , K 3 ). The expression for the mass operator is now 
simply 

M = E 12 + E 3 , M 3 = E 1 + E 2 . (49) 
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The use of light front variables enables one to separate the center of mass motion from 
the internal motion. The internal wave function \I/ is therefore a function of the relative 
momenta q and K. The function ^ is a product ^ = with $ = flavor, \ = spin, and 

<p = momentum distribution. The color wave function is antisymmetric. 

The angular momentum j can be expressed as a sum of orbital and spin contributions 



j = iV p x p + ^K M jS 



3 ' 



(50) 



3=1 



where 1Z m is a Melosh rotation acting on the quark spins Sj, which has the matrix repre 
sentation (for two particles) 

m + £M — ia ■ (n x q 



(X'\n M ^q ± ,m,M)\X) 



m + £M) 2 + q 2 ± 



(51) 



A' A 



with n = (0,0,1). The effects of the Melosh rotation are to significantly increase the 



computed charge radius |48 



The operator j commutes with the mass operator M; this is necessary and sufficient for 
Poincare-invariance of the bound state. In particular, j 2 \^, T>= 3/4|\l/, |> and j z \^, T> = 
l/2|\l/, |>. The angular momentum operator is in terms of relative coordinates given by 



j = iV K x K + TZ M (v, K ± , M 3 , M)ji2 + U M {1 - r], -K ± , m 3 , M)s 3 



ji2 = «V q x q + 1Z M {£, q±, m 1 , M 3 )s 1 + TZ M {1 - £, -q±, m 2 , M 3 )s 2 . 



(52) 



The orbital contribution does not contribute for the ground state baryon octet, so that 



Tin 



1 



ac - q R K L -aq L - cK L ' 



a? + K\\l c 2 + q\ \ cK R + aq R ac - q L K R 
I i ad + q R K L aq L - dK L 



a 2 + K\\j d 2 + q\ \ dK R - aq R ad + q L K 
b Kr 



(53) 



R , 



b 2 + K 2 \-K f 
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with 



a = M 3 + r]M, b = m 3 + (1 - rj)M , 
c = mi +£M 3 , d = m 2 + (1 - £)M 3 , 
= qi + H2 , QL = qi- m , (54) 
K R = K l + iK 2 , Kl — K\ — %K 2 . 

The momentum wave function can be chosen as a function of M to fulfill the requirements 
of spherical and permutation symmetry. The S'-state orbital function (f>(M) is approximated 
by either 

~ M 2 



0(M) = iVexp 



or ^ = WTWF^ (55) 



which depend on two free parameters, the constituent quark mass and the confinement scale 
parameter (3. The first function is the conventional choice used in spectroscopy, but it has 
a too strong falloff for large values of the four-momentum transfer. We use Schlumpf's 
parameters (3g=0.5Q GeV, (3=0.607 GeV, and the constituent quark mass, mj=0.267 GeV. 
The total wave function for the proton is given by 



p = —= (uudx X3 + udux X2 + duux X1 ^j 4> , 



V3 

(56) 



x f = — (ITT + TIT-2TT!), 



with 

1 

7 

x f = -L(2 lit - IT! - Tl!) • (57) 

The spin wave functions x X2 an d X X1 are the appropriate permutations of x X3 - The spin- wave 
function of the ith quark is given by 

T=7^ Q and l=n t . (58) 
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We now turn to the calculation of the proton form factors. The electromagnetic current 
matrix element can be written in terms of two form factors taking into account current and 
parity conservation: 



(N, \'p'\J»\N, \p) = u x ,(p') 



«a(p) 



(59) 



with momentum transfer Q 2 = — (p' — p) 2 and J M = g7 M g. For Q 2 = the form factors i*\ 
and _F 2 are respectively equal to the charge and the anomalous magnetic moment in units 
e and e/M^, and the magnetic moment is fi = Fi(0) + ^(O). The Sachs form factors are 
defined as 



O 2 

G z = F i-JT7T F 2> and G M = F 1 + F 2 
AM 2 N 



(60) 



and the charge radii of the nucleons are 
JF(Q 2 



-6- 



dQ 2 



and (r 2 E/M 



6 dG E/M (Q 



G E/M (0) dQ 2 



(61) 



The form factors can be expressed in terms of the + component of the current: 

Fi(Q 2 ) = ^(na\j + \na) , 



(62) 



Q±F 2 (Q 2 ) 



2M, 



N 



2P^ 



(iV,t|j + |iV,|) . 



The form factors are calculated from the diagrams of Fig|| using the "good" current J H 
so that no terms with qq pairs are involved. Schlumpf 's result is 



Fi(Q 2 



\E 3 E 12 M' 



(27T) 



1/2 



t (M')0(M) 



x 



i=l 



(63) 



Q±F 2 (Q 2 ) = -2M 



A? 



0T(M')0(M) 



i=l 
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with i = (uud) for the proton. Here the prime indicates the absorption of the momentum 
transfer as K' ± = K± + r/Qj_ and q' ± — q±. The factors F lu and F ld are the charges of the u 
and d quarks. We also consider Ga{Q 2 )- We take the hadronic axial- vector current to be 

A* = u'j^d . (64) 

The form factor of interest is given by 

2P + G a (Q 2 ) = (b',\ A + £?,T) • 



(65) 

The final step is to specify the operator r 2 for this relativistic model. We use 

r 2 = 3(p 2 + A 2 ) (66) 

— * — * 

where p and A are canonically conjugate to the momenta q and K: 

P=^^ A = -^(^ 1 + (l-Or 2 -f 3 ). (67) 

— * 

Note that p and A reduce to the usual three-body variables in the non-relativistic limit of 

VI. RESULTS 

It is worth while to begin by discussing whether or not point-like configurations occur in 
the relativistic models we employ. We do this by defining a quantity r 2 (Q 2 ) as 
2^2, (iV,T |r 2 J+|iV,t) 



r\Q 2 ) = 



(iv,t |J+|iv,t) 

I d3qd3K [flT^) E^(^'),xf|r 2 |xf,0(M)). (68) 



(2vr)6F 1 (g 2 ) J \E 3 E 12 M', 

This quantity should not be confused with the transverse size, b 2 (Q 2 ), used in studies of color 
transparency. The production of a point-like configuration at large momentum transfer by 
the electromagnetic current is signaled by the vanishing of the transverse size for large Q 2 . 
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Nevertheless, since we are interested here in the properties of a nucleon which is not necessar- 
ily moving at a large momentum with respect to the medium, it is r 2 (Q 2 ), as defined in (pBD, 
which enters our calculations. For comparison, we may define the quantity b 2 (Q 2 ) from (p5[) 
by including only a single transverse component of the operator r 2 . We see from Fig. |^ that 
both the power-law and the Gaussian wave functions display significant reductions of the 
quantity r 2 (Q 2 ) for increasing Q 2 , however only the power-law wave function is suppressed 
in the case of the transverse size, b 2 (Q 2 ). This behavior has been noted previously [^6|, and 
exemplifies the feature that power-law wave functions have point-like configurations in the 
transverse variables, but that such are absent in Gaussian wave functions. We shall use only 
the power law form, unless otherwise noted, because of its ability to reproduce data. 

Comparisons between the free form factors and the medium modified ones are given in 
Figs. ll^jQ- The quantity Ge is suppressed at low momentum transfer, but Gm is not. This 
is due to the relativistic nature of our model. Indeed, the magnetic moment is increased 
by about 5%. This small change is not enough to cause disagreement with existing nuclear 
phenomenology. At higher momentum transfers both form factors are inhibited by about 
10%. This is a significant effect, but not large enough to disagree with Jourdan's analysis 

0- 

This is more clearly seen by computing the nuclear response functions for the inclusive 
(e, e') cross section. (See for example @.) The excitation energy is u and the three 
momentum transfer is q so that 



"^L(q,^) + (^ + tan 2 |)i? T ( q ,^)] , (69) 



d 2 a 



0"M 



dQdE 

with the Mott cross-section a M = a 2 cos 2 (6 / 2) /4E 2 sin 2 (6/ 2). Here Q 2 = -g 2 = q 2 - u 2 
and 9 is the scattering angle. The longitudinal Rl and transverse Rt response functions are 
calculated in the relativistic Fermi gas approximation at a density p = 2k F /37r 2 , 

R L = -—Im(ZU p 00 + NU^) , (70) 

7Tp 

Rt = -±I m (ZU p 22 + NW 22 ) , (71) 
25 



for a target with Z protons and N neutrons. (Note, q is assumed to be along the 1 axis so 
the subscript 22 refers to a transverse direction.) Here the Fermi momentum is taken to be 
kp = 260 MeV, which is appropriate for 56 Fe. 
The polarization II 

d A p 



(27T) 



r Tr[G(p + g)rG(p)rj 



is calculated with the nucleon Greens function G(p) 



G(p) 



2E n 



Q(k F -\ P \) , e(\p\-k F ) 



(72) 



(73) 



_p - E p - ie po - E p + it 
and the electromagnetic vertex T, for i = p (proton) or n (neutron). We assume the elec 
tromagnetic vertex has the form, 



^-^ + ^2 2M 



(74) 



even for off-shell nucleons in the medium. We show results comparing using Fi computed in 



free space and in the medium in Figs. 12 and 13 



It is immediately evident from Fig.[T3 that modifications to the magnetic form factor 
Gm have no effect on the longitudinal response, while the modifications to the electric form 
factor Ge lead to a suppression. This effect is due to the larger charge radius in the medium. 



The situation is quite different in the case of the transverse response shown in Fig.[13|. There 
it is seen that modifications to the electric form factor have no effect, while modifications to 
the magnetic form factor lead to only a minor suppression. In this case the increase in the 
magnetic moment, as shown in Fig.|Tl~|, tends to cancel the effect of the increased radius. 



We also study the effects of using the form factors of Figs. |10| and [11] in computing the 
(e, e'p) cross sections for finite nuclei. In this case, absorption effects emphasize the role of 
the nucleon surface and the influence of color neutrality is about 60% smaller than shown 



in Figs. 12 and 13 



Our results for the medium modifications on the form factor Ga are shown in FigJTJ. 
Once again there is about a 10% reduction. This effect could be observed in neutrino-nucleus 
scattering, in the (p,n) reaction, or parity violating electron scattering ||65|| . 
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We now turn to computing the valence structure functions in the medium and in free 
space. One advantage of the light front formalism is that the wave function is closely related 
to the valence structure function F% al : 

FrKxBj) = J d 3 qd 3 K^5 (x Bj - |M . (75) 

Thus the operator O of Eq. (|l^) is the delta function which sets the plus-momentum of 
the quark equal to xsjP + - This is meant to be taken at some momentum scale ~ 1 GeV. 
We may compute this quantity for the free and medium modified wave function. There 
are a host of other effects of the medium including Fermi motion, shadowing, pions in the 



medium, nuclear correlations, six quark bags, etc. (see the reviews |T(J.) Our concern here 
is to assess the effects of color neutrality. So we do not include these effects and do not 
compare our results with data. 



The free version is shown in Fig.|15|, where the expected shape is obtained. The results 
for the ratio of medium modified to free structure functions using both the power law and 
Gaussian forms are shown in Fig.jTBI. This shows that suppression of point like configurations 



does indeed lead to suppression of the valence structure function at large values of x B j- The 
normalization of the wave function ensures that the integral / dxBjF 2 (xBj) is unity whether 
the free or modified wave function is used. Thus one expects regions in which the ratio is 
bigger and smaller than unity. However, the result that there is suppression at large XBj is 
not a trivial consequence of kinematics and normalization as is seen by comparing with the 
Gaussian form. The Gaussian wave function does not display as much suppression at large 
Xsj, which indicates the relative importance of large momentum components, or point-like 
configurations, in the power law versus Gaussian wave functions. The suppression occuring 
at large XBj is consistent with the relevant features of the data and therefore provides 
evidence for the existence of point like configurations. 
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VII. SUMMARY AND DISCUSSION 



We have used the ideas of color neutrality to motivate a functional form of the central 
shell model potential: 



V(r,R) = V oP (R)—^-, (76) 



< r 2 > 



where < r 2 >=< 0ol r2 |0o > is the nucleonic expectation value of the operator r 2 , Vq ~ 
— 50MeV and p(R) is the nuclear density normalized so that p(R = 0) = 1. The concern 
here is with fluctuations, so that we rewrite the central shell model potential as 



2 2 

V(r,R) = V p(R)+V p(R) r "y^ , (77) 

and treat the second term as a perturbation. The effects of this perturbation can be evalu- 
ated, using Eq.flTTD (or using the closure approximation of Eq. flT^) ) for any nuclear process. 

Solving a set of toy models indicates that a perturbative treatment is valid and further- 
more that a closure approximation may be used to avoid computing the complete spectrum 
of baryonic wavefunctions usually necessary in perturbation theory. The toy-model results 
are consistent with the notion that point like configurations can be suppressed in the nuclear 
medium. 

Light front quantum mechanics, along with a specific model of the nucleon wave function 
34H36|, is next employed to compute the influence of medium effects at relatively high 



momentum transfer. We find that at low values of Q 2 < 1 GeV 2 the electric form factor is 
suppressed and displays an increased charge radius, but that while the magnetic radius is 
also increased so is the magnetic moment. This leads to the result that the (e, e') transverse 



response, shown in Fig.|13], is largely unaffected by the medium in this context. This behavior 
of the transverse response has previously been interpreted as signaling no change in the 
magnetic radius -contrary to the result obtained here. At higher values of Q 2 both form 
factors are suppressed in the medium as is Ga- These results are in accord with ideas about 
the suppression of the longitudinal response, but are not inconsistent with the analysis of 
Jourdan |13 |. 
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The results for the medium modifications of the F2 structure function show a suppression 
at large values of XBj for both the power-law and Gaussian wave functions. However, the 
suppression is greater in the case of the power law, indicating the more pronounced role of 
high momentum components or point-like configurations there. This is consistent with the 
works of and p3| , pT| and more importantly, provides evidence for the existence of point 



like configurations. 
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FIGURES 

FIG. 1. Possible configurations of the proton wave function. 

FIG. 2. The matrix element of Eq. (12) (j„ commutes with h™ 1 .) 

FIG. 3. The influence of the attractive Coulomb potential in Eq.(^) on the wave function is 
illustrated. 

FIG. 4. The influence of the medium on the wave function is illustrated as described in the text. 

FIG. 5. The medium-modified exact wave function and that calculated in the closure approxi- 
mation are plotted as a function of the relative coordinate, r. 

FIG. 6. The exact medium-modified, free, and the closure-approximated form factors are plot- 
ted as a function of the square of the momentum transfer. 

FIG. 7. Medium modified and free form factors calculated from the analytic expression given 
in Eq.(f4l|) are plotted as a function of the square of the momentum transfer. 

FIG. 8. The absorption of momentum by the valence quarks is illustrated. 

FIG. 9. The quantities r 2 (Q 2 ) and b 2 (Q 2 ) for power-law and Gaussian wave functions. 

FIG. 10. The free and medium modified electric form factors vs. the square of the 
four-momentum transfer, Q 2 . 

FIG. 11. The free and medium modified magnetic form factors vs. Q 2 . 

FIG. 12. The longitudinal response vs the energy transfer, u>. 
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FIG. 13. The transverse response vs. cj. 



FIG. 14. The free and medium-modified axial vector form factors vs Q 2 . 

FIG. 15. Proton structure functions are plotted as a function of the scaling variable xbj for 
power-law and Gaussian wave functions. 

FIG. 16. The ratio of the medium-modified to free structure function is plotted as a function 
of the scaling variable for the power-law and Gaussian wave functions. 



37 



